Upper critical field for anisotropic superconductivity. A tight— binding approach. 
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We study the problem of the upper critical field (H C 2) 
for tight-binding electrons in a two-dimensional lattice. The 
external magnetic field is introduced into the model Hamilto- 
nian both via the Peierls substitution and the Zeeman term. 
Carrying out calculations for finite systems we analyze the 
influence of the external field in the commensurable and in- 
commensurable case on an equal footing. The upper critical 
field has been discussed for intrasite as well as anisotropic 
intersite pairing that, in the absence of magnetic field, has a 
d x 2_ y 2 symmetry. A comparison of H C 2 determined for differ- 
ent symmetries shows that the on-site pairing is more affected 
by the external field i.e., the critical temperature for the on- 
site pairing decreases with the increase of the magnetic field 
faster than in the anisotropic case. Moreover, we have shown 
that the tight-binding form of the Bloch energy can lead to 
the upward curvature of H C 2 , provided that the Fermi level is 
close enough to the van Hove singularity. 
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I. INTRODUCTION 

One of many striking properties of high -temperature 
superconductors is related to the field-induced transition 
from superconducting to normal state. Magnetic proper- 
ties of high-T c compounds give rise to both quantitative 
and qualitative differences with respect to the conven- 
tional superconductors. The systems under considera- 
tion are characterized by extremely high values of the 
upper critical and its unusual temperature dependence. 
For optimally doped samples experimental investigation 
of the critical field is limited only to temperatures close 
to T c jij, whereas at lower temperatures the magnitude 
of H C 2 is far beyond the reach of laboratory magnetic 
fields. The measurements carried out in a wide range 
of temperature for underdoped superconductors clearly 
indicate the positive curvature of H c i{T) even at gen- 
uinely low temperatures pJ-pUj. Theoretical approaches do 
not provide a unique, complete description of these phe- 
nomena. The most of unconventional properties of high- 
temperature superconductors, like narrow quasiparticle 
bands, lifetime effects of states close to the Fermi level 
and linear temperature dependence of the normal-state 
resistivity are usually attributed to strong Coulomb cor- 
relations. However, upward curvature of the upper criti- 
cal field is observed also in overdoped compounds, where 
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the temperature dependence of resistivity changes gradu- 
ally from linear to quadratic behavior . This feature 
suggests that the positive curvature of H C 2 (T) could orig- 
inate from e.g., symmetry of the superconducting order 
parameter or details of the density of states and may be 
explained without a sophisticated treatment of the most 
difficult problem that is related to the presence of strong 
electronic correlations. 

It is believed that the symmetry of the superconduct- 
ing state can be close related to the pairing mechanism. 
There is a lot of node-sensitive experiments, based on the 
angle resolved photoemission spectroscopy [Q, London 
penetration depth [|| , NMR || and quasiparticle tunnel- 
ing pp|| , which indicate that the energy gap is strongly 
anisotropic and vanishes in particular directions in the 
Brillouin zone. Moreover, the phase-sensitive super- 
conducting interference device experiments pd| ] demon- 
strated the sign change of the order parameter between 
the x and y directions. Generally, these results are consis- 
tent with the d x 2_ y 2 pairing scenario. On the other hand 
there are experimental indications, which had questioned 
the pure d x 2_ y 2 symmetry of the energy gap and suggest 
mixed pairing symmetry with a dominant d-wave com- 
ponent (e.g., d ± s or d ± is) ]f^-|il| . 

The measurement of the upper critical field can give 
insight into the microscopic parameters of a relevant 
model. For example, the coherence length £ is usually 
derived indirectly from the expression H C 2(0) = </>o/27r£ 2 
[PH , where H c2 (0) is the upper critical field determined 
at T = 0, and 0o is the magnetic flux quantum. The 
theoretical investigation of the upper critical field for 
different pairing symmetries is predominantly based on 
the Ginzburg-Landau (GL) [[l7| theory or the Lawrence- 
Doniach |18| approach in case of layered superconduc- 
tors. With the help of linearized GL equations Won and 
Maki [|l9| have shown that H C 2 in a model with repulsive 
on-site interaction depends linearly on temperature near 
T c and saturates at T — > 0. They have not found any 
sign of the upward behavior. There are also calculations 
for H C 2 in systems with mixed symmetries, especially for 
superconductors in which the dominant ci-wave order pa- 
rameter coexists with a subdominant s-wave component. 
However, in the most of these approaches H C 2{T) exhibits 
negative curvature. On the other hand, results obtained 
in Ref. [^o) suggest that the upward curvature of the 
critical field could be a characteristic feature of a d-wave 
superconductor. The positive curvature of H c2 (T) can 
also originate from the presence of magnetic impurities 

HP- 

A separate problem, that is usually neglected in the 
above approaches, is the influence of the periodic lattice 
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potential on the upper critical field p3[ . Application of 
magnetic field to the two-dimensional (2D) electron sys- 
tem in a tight-binding approximation leads to a fractal 
energy spectrum known as Hofstadter's butterfly, where 
even very small changes in magnetic field can result in 
a drastic changes of the spectrum |24| |2(| . In this paper 
we investigate the upper critical field for electrons de- 
scribed by the two-dimensional tight-binding model with 
intra- and intersite pairing. We show that anisotropic su- 
perconductivity is less affected by the external magnetic 
field then the isotropic one. We also demonstrate that 
the lattice effects can give rise to important corrections 
with respect to Helfand-Werthamer |^7],|lJ solution of 
the Gor'kov equations pq] . This effect is of particular 
importance in the vicinity of the van Hove singularity 
and at low temperatures. 

II. GAP EQUATION CLOSE TO H C2 

We consider a two-dimensional square lattice im- 
mersed in a uniform, perpendicular, magnetic field. The 
BCS-type Hamiltonian is of the form 

H = H^in +Hy - /X ^ c^Cjcr 

-9HbH z ( c ! T c *T - 4i c *i) > (!) 

i 

where c\ a (ci a ) creates (annihilates) an electron with spin 
a on site i. The chemical potential /j, is introduced in or- 
der to control the doping level. The last term in the above 
Hamiltonian describes the paramagnetic Pauli coupling 
to the external field. Here, g stands for the gyromag- 
netic ratio, [1b is the Bohr magneton and H z is the z— 
component of the external field. The first (-ffkin) an d 
the second (Hy) term in the Hamiltonian represents the 
kinetic energy and the pairing interaction, respectively. 
Within the tight-binding approach 

#kin = J2 M A KW- (2) 

The electrons are gauge-invariantly coupled with local 
U(l) gauge field by a phase-factor in the kinetic-energy 
hopping term. According to the Peierls substitution |2jJ 
in the presence of magnetic field the original hopping 
integral between sites i and j, tij acquires an additional 
factor 

tij (A) = Uj exp(^-J^ A ■ dlj , (3) 

In the case of the on-site pairing, which leads to isotropic 
order parameter, the BCS-type interaction takes on the 
form 



Here, we have introduced local superconducting order 
parameter, A, = (c^Qf), which in the presence of the 
magnetic field can change from site to site . We also 
consider anisotropic superconductivity with the intersite 
pairing interaction given by 

H V = -VJ2 (4A A « + A tt) • (5) 

<ij> 

For the sake of simplicity we restrict our considerations 
only to the nearest-neighbor coupling with the singlet 
order parameter Ay = (cix c iT — c «T c ii)- 

We start with the discussion of the normal state prop- 
erties. Similarly to Ref. |25| ] we make use of an unitary 
transformation U that diagonalizes the kinetic part of the 
Hamiltonian 

U^H kin U = tt kin . (6) 

This transformation defines a new set of fermionic oper- 
ators a na = J2i U^iCicr, m which the Hamiltonian in the 
normal state takes on the diagonal form 

U = ^ i E n - M - crgfi B H z )al l<T a ncr . (7) 

li a 

In the absence of the magnetic field U represents trans- 
formation from the Wannier to the Bloch representation. 
For finite magnetic field and general gauge the quantum 
number n enumerates eigenstates, although does not rep- 
resent a reciprocal lattice vector. In order to simplify 
further discussion we restrict our considerations only to 
the nearest neighbor hopping with t^j\ = —t. We also 
assume the type-II limit of superconductors where the 
magnetic field can be regarded as a spatially uniform ob- 
ject. Choosing the Landau gauge A = H z (0, x, 0) the 
hopping integral depends explicitly only on x and the 
momentum in y direction p y remains a good quantum 
number. Due to the plane-wave behavior in y direction 
the unitary matrix U takes on the form 

U l(P ^ Py) = U (Xty){Px , Py) = TV" 1 / 4 e i Pyy a g(p x , Py ,x) , (8) 

where (ax, ay) is the position of the i-th site and 
(p x ,p y ) represents n-th eigenstate of the Hamiltonian 
(7). Straightforward calculations ||(| show that the x- 
dependent part of the wave function g (p x ,p y , x) fulfills a 
one-dimensional difference equation 

g (p x ,p y , x + 1) + 2 cos (hx - p y a) g (p x ,P v ,x) 
+g(Px,P y ,x- 1) = t^ 1 E(p x ,py)g(p x ,p y ,x) , (9) 

where we have introduced a reduced dimensionless mag- 
netic field h = ea 2 H z /(hc). This quantity can be ex- 
pressed with the help of magnetic flux <fi through lat- 
tice cell and flux quantum (h = 2it(j)/An). Equation 
(9) is known as the Harper equation [ p4[ and has ex- 
tensively been studied . The Harper equation, de- 
rived here within a tight-binding approximation, can be 
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also obtained in a case of weak perturbation of a Landau- 
quantized two-dimensional electron system |25| , ^] . 

Now, let us take into account the pairing potential Hy . 
In order to investigate the transition from the supercon- 
ducting to the normal state we make use of equation of 
motion for the anomalous Green function. In the case of 
the isotropic on-site pairing one obtains 

[v - E(p x , Py )+ fi + gfi B H z ] {{a (PxtPy)1 | a(^ jfey ) i )) 
= ~ V E ^%*)^)<<^);l«(S.,*,)l>>- 



where 



(10) 



As far as we are close to the phase transition we make 
use of a linearized gap equation i.e., we calculate the 
propagator (( a |^, k , ^ I a (^. k ^)) m the normal state. 

Similarly to the standard BCS theory, such approach al- 
lows one to determine the critical temperature or, in our 
case, the upper critical field. However, it is irrelevant for 
calculations below T c . 

The choice of the Landau gauge implies that the 
isotropic order parameter does not depend on y: Aj = 
^(x,y) = Then, the linearized gap equation reads 



A = MA, 
where A = (Ai, A 2 , A3, ...) and 



(11) 



M(x,x') = -j= 9{Px,P v ,x) g(k x ,-p y ,x) 

x 9 ( h , ~Py ,x') g(p x ,p y ,x') x(Px,P y ;k x ,-p v )- 

(12) 

In the presence of the magnetic field the Cooper pair 
susceptibility is given by 



x{Px, Py] k x , ky) 



tanh 



E(p x ,p y ) - p - 9HbH z 



tanh 



2k B T 

E {k Xl k y ) - p, + g[i B E z 



2k B T 

x [2(E{p x ,p y )+E(k x ,ky) -2/i)] _1 . 

(13) 

In the case of the nearest-neighbor pairing we obtain 
the gap equation analogous to Eq. (11). Similarly to 
the isotropic pairing Ajj does not depend explicitly on y. 
However, there are two types of order parameter at each 

(x) 

site: A x when sites i and j lay along the x axis, and 
Ax when sites i and j lay along the y axis. Close to 
the upper critical field the gap equation for anisotropic 
superconductivity can be written in a matrix form 



/AW 
^ A^) 



(14) 



V 



MS a ^{x,x') = -j= x(px,p y ;k x ,-p y ) 

Px iPy j^x 

xA (q) (p Xl k Xl p yi x) A {(3) (p x ,k x ,p y ,x') , 

(15) 

and 

A {x) (p x ,k x ,p y ,x) = g(p x ,py,x)g(k x ,-p y ,x + 1) 
+ g(px,p y ,x+ l)g(k x ,-p y ,x), 

(16) 

A {y) (p xi k x ,p y ,x) = 2cos(p y a) 

x g(Px,P y , x)g(k x , -p y , x). (17) 

Equations (11) and (14) constitute a system of linear 
equations for the order parameters and the condition for 
existence of a non-zero solution can be written as 



det (M - I) = Q 
in the case of isotropic pairing, and 



det 



J^(x,x) _ j J^(x,y) 

M^ V,X ^ M^ y ' y "> - 1 



(18) 



(19) 



for anisotropic superconductivity, where / is the unit ma- 
trix. These equations allow one to obtain the magnitude 
of the upper critical field perpendicular to the plane. For 
the two-dimensional square lattice the size of matrices 
which enter Eqs. (18) and (19) is proportional to the 
square root of the number of the lattice sites. Analyti- 
cal solutions of the Harper equation (9) are known only 
in a few cases of commensurable field pl| (in our nota- 
tion h — 2np/q, where p and q are relative prime inte- 
gers), which correspond to unphysically high magnetic 
field. Therefore, in order to investigate H C 2 we restrict 
our considerations to a finite lattice, for which we are 
able to analyze numerically the commensurable and in- 
commensurable magnetic field on an equal footing. 



III. DISCUSSION OF RESULTS 

We consider square M X M cluster with periodic 
boundary conditions (be) along the y axis. As the Lan- 
dau gauge breaks the translation invariance along x axis 
we use fixed be in this direction. An additional advan- 
tage originating from such a mixed be is the absence of 
the unphysical degeneracy of states at the Fermi level, 
which occurs for the half-filled band in cluster calcula- 
tions with fixed or periodic be taken in both directions 
J33|. In order to estimate the finite size effects we have 
carried out numerical calculations for clusters of differ- 
ent sizes. We have found that in the case of the isotropic 
pairing and small concentration of holes (8 < 0.2) there 



3 



are no significant differences between results obtained on 
150 x 150 and 200 x 200 clusters. For anisotropic pairing 
already 120 x 120 clusters give convergent results. 




0.001 



0.000 



0.000 



0.010 



0.020 



0.030 



k B T/t 



FIG. 1. Temperature dependence of the reduced upper crit- 
ical field for isotropic pairing and different occupation num- 
bers n. The cross, circle and square marks indicate results ob- 
tained on 150 x 150 cluster, whereas the solid lines correspond 
to 200 x 200 cluster. The arrows show the superconducting 
transition temperature for an infinite system calculated from 
the BCS gap equation in the absence of magnetic field. V = t 
has been assumed. 



order parameter i.e., for isotropic (Fig. 1.) as well as 
anisotropic pairing (Fig. 2.), the slope of H c2 (T) strongly 
decreases with increasing doping. Note that our cluster 
results exactly reproduce the BCS transition tempera- 
ture when the magnetic field tends to zero. In the case 
of intersite pairing the arrows indicate BCS solutions for 
d x 2_ y 2 superconductivity. However, the external mag- 
netic field affects the relative phases of the order param- 
eter in the x and y directions, which can change from 
site to site. Therefore, it is impossible to determine glob- 
ally the type of the symmetry of the energy gap in the 
presence of magnetic field. 

Contrary to the conclusion presented in Ref. [^0| , our 
results (Figs. 1. and 2.) do not indicate that the upward 
curvature of H C 2 (T) can emerge as a direct consequence 
of the symmetry of superconducting state. However, the 
anisotropy of the order parameter can significantly influ- 
ence the magnitude of the upper critical field. In order 
to investigate this relationship we have directly compared 
results obtained for on- and intersite pairing for isotropic 
and anisotropic superconductivity. We have chosen the 
magnitudes of the pairing potentials V, which, in the ab- 
sence of magnetic field, lead to the same superconduct- 
ing transition temperatures for isotropic and anisotropic 
superconductivity. Fig 3. shows the temperature depen- 
dence of the upper critical field obtained for the half-filled 
case. One can see that the anisotropic superconductivity 
is less affected by the external field than the isotropic 
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FIG. 2. The same as in Fig. 1. but for anisotropic pairing. 
The cross, circle and square marks indicate results obtained 
on 120 x 120 cluster, whereas the solid lines correspond to 
150 x 150 cluster. The arrows show the d-wave supercon- 
ducting transition temperature for an infinite system calcu- 
lated from the BCS gap equation in the absence of magnetic 
field. Here, V = 0.3 1 has been assumed. 

Figures 1. and 2. show the reduced critical field, 
h C 2 = ea 2 H C 2/(hc), for different concentrations of holes. 
Independently on the symmetry of the superconducting 
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FIG. 3. The temperature dependence of the upper critical 
field evaluated for the half-filled case, n = 1. The circle and 
diamond symbols denote results obtained for intersite pairing 
with V = 0.244t and on-site pairing with V = t, respectively. 
The symbols and continuous lines correspond to the same 
sizes of clusters as in Figs. 1 and 2. The inset shows the upper 
critical field obtained for the occupation number n = 0.8. 
Here, we have taken the intersite pairing potential V = Q.3t. 

An important observation is that this result depends nei- 
ther on the magnitude of the pairing potential nor on the 
concentration of holes (see the inset in Fig. 3). Thcrc- 
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fore, it can be considered as a characteristic feature of 
the two-dimensional lattice gas. 

In the absence of magnetic field there is a van Hove 
singularity in the middle of the band. Although, the ex- 
ternal field results in a splitting of the Bloch band into 
a huge number of subbands, the presence of the origi- 
nal van Hove singularity is reflected in the Hofstadter 
spectrum |2q| . In contradistinction to the structure of 
Landau levels, the Hofstadter spectrum does not consist 
of uniformly distributed energy levels. In particular, the 
average distance between the energy levels close to the 
Fermi energy achieves its minimum when the chemical 
potential is in the middle of the Bloch band. It can be 
considered as a remnant of the original van Hove singu- 
larity. 
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FIG. 4. The upper critical field for the on-site pair- 
ing (dashed lines) fitted to the results obtained for the 
two-dimensional version of the Helfand-Werthamer approach 
to the Gor'kov equations (continuous lines). We have chosen 
the pairing potential V = t, for different values of the occu- 
pation number n, indicated in the figures. The inset shows 
results obtained for V = 0.7t and n = 1. 

The question which arises concerns the impact of this 
feature on the upper critical field. In order to analyze 



this problem we have fitted H C 2 (T) obtained for isotropic 
superconductivity to the results obtained for the two- 
dimensional version of the Helfand-Werthamer ap- 
proach to the Gor'kov equations. Fig. 4 shows the 
numerical results. Away from the half-filled case the 
qualitative temperature dependence of the upper critical 
field can be very well approximated by the solution of 
the Gor'kov equations. It suggests, that the complicated 
Hofstadter spectrum does not influence the temperature 
dependence of the critical field, provided that the Fermi 
level is far enough from the original van Hove singular- 
ity. However, in the vicinity of the van Hove singularity 
the second derivative of H C 2(T) is significantly enhanced, 
when compared to the results obtained from the Gor'kov 
equations. It is of particular importance for small values 
of the pairing potential, when the system remains in su- 
perconducting state only at relatively low temperatures 
and the Cooper-pair susceptibility is strongly peaked at 
the Fermi level. Then, the curvature of H C 2 (T) can grad- 
ually change from negative to positive, as depicted in the 
inset in Fig. 4. This effect takes place for isotropic as 
well as for anisotropic pairing. Similar results have been 
reported in Ref. [B4J. 



IV. CONCLUDING REMARKS 

We have investigated the temperature dependence of 
the upper critical field for the two-dimensional lattice 
gas. With the help of unitary transformation we have 
obtained a diagonal form of the Hamiltonian in the nor- 
mal state and derived gap equations both for isotropic 
and anisotropic superconductivity. We have discussed 
influence of the symmetry of the superconducting state 
and the van Hove singularity on the upper critical field. 
Our results clearly indicate that the symmetry of the su- 
perconducting order parameter itself can not lead to up- 
ward curvature of H C 2(T). However, quite pronounced 
tendency can be observed for the half-filled case, when 
the Fermi energy is close to the original van Hove singu- 
larity. In the absence of the external field this singularity 
occurs in the middle of the band. The enhancement of 
curvature of H&iT) takes place for isotropic as well as 
anisotropic superconductivity and is of particular impor- 
tance for small values of the pairing potential. Then, the 
curvature can gradually change from negative to posi- 
tive. This effect smears out for larger doping where the 
temperature dependence of the upper critical field can 
be rendered very well when solving the Gor'kov equa- 
tions. We have found that in the case of anisotropic 
pairing the upper critical field exceeds the critical field 
obtained for isotropic superconductivity. It takes place 
for small doping [6 < 0.2) and arbitrary magnitude of 
the pairing potential. These results suggest that in the 
two-dimensional lattice gas anisotropic superconductiv- 
ity is less affected by the external field than the isotropic 
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The proposed method allows one to derive the gap 
equation in the same way as the standard BCS approach. 
The only differences are related to the fact that the diago- 
nal form of the normal-state Hamiltonian is obtained nu- 
merically and the superconducting order parameter can 
be a site-dependent quantity. The similarity between our 
method and the BCS approach allows for straightforward 
incorporation of the local Coulomb repulsion within any 
standard approximation. Here, one may expect destruc- 
tive influence of correlations, in particular in the isotropic 
channel. This originates from the fact that local repul- 
sion always acts to the detriment of the formation of 
local Cooper pairs. The impact of Coulomb, Hubbard- 
like, correlations on anisotropic superconductivity seems 
to depend on the approximation scheme. This problem 
is under our current investigation. 
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